The Field Theory of the q -> 4+ Potts Model 



O 



G. Delfino a and John Cardy 5 ' c 

a Laboratoire de Physique Theorique et Hautes Energies 

Universite Pierre et Marie Curie, Tour 16 l er etage, 4 place Jussieu 

75252 Paris cedex 05, France 

^ ■ b Theoretical Physics, University of Oxford 

^ ' 1 Keble Road, Oxford 0X1 3NP, United Kingdom 

PL, " 

in 



All Souls College, Oxford 



> 

^—i | Abstract 
(N 

The g-state Potts model in two dimensions exhibits a first-order transition for q > 4. 
As q — > 4 + the correlation length at this transition diverges. We argue that this 
limit defines a massive integrable quantum field theory whose lowest excitations are 
kinks connecting 4 + 1 degenerate ground states. We construct the S-matrix of this 
theory and the two-particle form factors, and hence estimate a number of universal 
amplitude ratios. These are in very good agreement with the results of extrapolated 
series in q^ 1 ^ 2 as well as Monte Carlo results for q = 5. 
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The g-state Potts model, defined by the lattice Hamiltonian [1| 

H = -J '^2 S s (x),s(y) > C 1 ) 

where the spin variable s(x) assumes q different values (colours), continues to be of fun- 
damental importance in the description of a large variety of critical phenomena, ranging 
from ferromagnetism to percolation and adsorbed monolayers [[| . It was shown by Baxter 
that in two dimensions the ferromagnetic model undergoes a phase transition which is 
continuous for q < 4 and first order otherwise ||. The Coulomb gas @] and conformal 
field theory (CFT) || provided later a complete description of the second order phase 
transition line, which turned out to correspond to CFT's with central charge c < 1, the 
value c = 1 corresponding to the end point q = 4. More recently, integrable field theory 
has lead to new results for the scaling limit of the off-critical model for q < 4 || [7], 

Concerning the first order transition for q > 4, several exact lattice results - internal 
energy [|J, magnetisation @, correlation length |I| - are known, but progress through 



field theoretic methods is generally prevented by the absence of a scaling limit. When q 
approaches 4 from above, however, the correlation length at T c diverges and a continuum 
description in terms of a massive quantum field theory should be possible. The identifi- 
cation and solution of the quantum field theory describing the limit q — > 4 + at T c , as well 
as the determination of the universal critical quantities, are the subject of this note. 



The correlation length at criticality of the lattice model is known to behave as |10 

f ~ ae^ 1 ^ (2) 

as q — > 4 + , where a is proportional to the lattice spacing, so that the continuum limit 
corresponds to taking the limits a — > 0, q — > 4 + in such a way that £ remains finite. The 
presence of an essential singularity rather than a power law divergence in Eq. (Q) is char- 
acteristic of a perturbing operator which is only marginally relevant (scaling dimension 
2) - except that in the original description of the Potts model q is a parameter which 
should be invariant along RG trajectories and therefore such an interpretation needs to 
be treated with some care. In fact we know that at q = 4 there is such a marginal field 
■0: it was shown by Nienhuis et al. [\ l| to correspond to the fugacity for vacancies in the 
lattice model. When ip < the transition is second order, with logarithmic modifications 
to scaling arising from the marginal irrelevance of ip, but when ip > 0, the transition is 
first order with a correlation length at the transition which diverges as ip —>■ + like 

£ ~ ae 1 ^ (3) 
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where b is a constant. 

The main point is that the scaling limits in which a — > with £ fixed are identical in 
these two cases. This may be seen, for example, from the general structure of the RG 
equations near q = 4 and ip = 0. Based on the analysis of Nienhuis et al. fli~ll , Cardy, 
Nauenberg and Scalapino |12| argued that these have, to lowest order, the general form 

dip/dl = bip 2 + a(g-4) + 0(^ 3 , (g-4)^,(g-4) 2 ) (4) 
dt/dl = (y t + ctij)t + 0(i} 2 t,(q-4)t,t 2 ) (5) 
dh/dl = (y h + c h i))h + 0{i) 2 h,(q-A)h,th) (6) 

where t is the deviation from the critical temperature and h is a symmetry-breaking field, 
a, b, Ct and are all constants, certain combinations of which are universal [|l2j. An 
important feature of (|5],|6|) is that, to lowest nontrivial order, they do not involve q. This 
is because, as may be seen from the first equation, q — 4 is effectively 0(ip 2 ). Integrating 

up to a value I such that if) (I) = 0(1) then gives results for the correlation length 
£ ~ e l in agreement with (§J^) in the two cases q > 4 and q = 4, ip > 0. The various 
thermodynamic quantities are then found by integrating the other equations up to / ~ ln£. 
To the order stated, the results will be identical in the two cases, when expressed in terms 
of £. Therefore the scaling limits are identical. 

Another way of understanding this is through the mapping of the lattice Potts model to 
a height model and thence to a Coulomb gas or sine-Gordon theory Q. From this point of 
view, q is merely a parameter identified with a certain function of the coupling constant 
conventionally called (3. q = 4 corresponds to f3 2 = 8tt at which point the operator 
corresponding to ip becomes marginal. Within the standard RG picture of the sine- 
Gordon model, both (3 and ip have non-trivial marginal flows. However the scaling limit, 
corresponding to the massive sine-Gordon theory at (3 2 = 8ir, is unique, and therefore 
describes both the cases q — * 4 + and ip — * + at q = 4. 

Having made this observation, we may take over the results of Ref . in which it was 
pointed out that the scaling limit of the massive q = 4 theory is integrable, and in which 
the scattering theory and form factors were determined. 

In our case, the construction of the scattering theory goes as follows. Along the first 
order phase transition line, q ordered ground states are degenerate with the disordered 
ground state. The field theory describing the scaling limit q — > 4 + has 4 ordered vacuum 
states Qi, i = 1, ... ,4. Invariance under colour permutations implies that, in the order 
parameter space, they lie at the vertices of a tetrahedron having the disordered vacuum 
Qq at its center. The elementary excitations of the scattering theory are stable kinks 
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Koi(9), Kio(6) interpolating between the center of the tetrahedron and the i-th vertex, 
and vice versa. We denote by 9 the rapidity variable parameterising the on-shell momenta 
as (jAp 1 ) = (m cosh 9, m sinh 9). The mass of the kinks m ~ C," 1 measures the deviation 
from the conformal point q — 4. The space-time trajectory of a kink on the plane 
draws a domain wall separating a coloured phase from the disordered one. The space 
of asymptotic states is made of multi-kink sequences in which adjacent vacuum indices 
belonging to different kinks have to coincide. For example, up to possible bound states, 
the lightest excitation interpolating between two ordered vacua is Kj a (8i)K j($2) . 

The factorisation of multi-kink processes reduces the scattering problem to the deter- 
mination of the two-kink amplitudes. Colour permutation symmetry allows only for the 
four elementary processes depicted in Fig. 1. The four amplitudes can be determined as 
a solution of the requirements of unitarity (crucially simplified by the absence of particle 
production), crossing and factorisation. The scattering amplitudes of Fig. 1 are given in 
[ PI and read 

P'^ 9 2?7T — 9 

Mo) = V ^Sr > ^ 

l m — o 

Ai(9) = ^—-?—S (9), (8) 

L ITX — U 

B Q (9) = e^ i ^S (9), (9) 
m — 9 

B 1 (9) = e** S (9) , (10) 
where 9 is the rapidity difference of the two kinks, 7 = - In 2, and 

( n\ F (I + 2§f) F (~2fe) (.[<*> dx e-f . x9\ 

S o( ) = ^rr > \ - / a v = - exp \ i I — —-jsm— >. (11) 



r(i-^)r(^) Wo x coshf 



The absence of poles in the physical strip Im# e (0, m) ensures that there are no 
bound states and that the four amplitudes above completely determine the scattering the- 
ory. This S-matrix shares evident analytic similarities with that of the S'C/(2)-invariant 
Thirring model. As a matter of fact, the latter is a realisation of the same perturbed CFT 
on a different particle basis (an SU(2) doublet rather than our kinks). The possibility for 
a single perturbed CFT to be invariant under different symmetry groups and to describe 



different universality classes is discussed, for example, in Ref. [14 



Making contact with the thermodynamics requires the computation of correlation 
functions. In our 5*-matrix framework these are obtained as spectral series summing over 
all multi-kink intermediate states. Neglecting terms of order e~ im ^ in this large distance 
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expansion, we will approximate the (connected) two-point correlator of a scalar operator 
$(x) as 

(fi |$(x)$(o)|Q )^E/ P-^tSVi-ettfe-M*, (12) 

i=l Jd 1 >e 2 27T Z7T 

in the disordered phase, and 

$k\*{x)m\su) * L f^l^i-^l 2 ^ 2 , (is) 

in the i-th ordered phase. Here E 2 = m(cosh.6i + cosh 8 2 ) is the energy of the two-kink 
asymptotic state and we introduced the two-kink form factors 

- 9 2 ) = (l} o |<KO)|i^(0i)^o(02)> , (14) 
F$(0i - 6 2 ) = (^|$(O)|if40i)K o *(02)> • (15) 

The operators of interest for us are the spin (Tj(x) = 5 a ( x )j — 1/?, and the energy 
£ ( x ) — I2yfis(x),s(y), whose scaling dimensions around the q = 4 fixed point are X a = 1/8 
and X £ = 1/2 @|. 

Some consequences for the physics of the coexisting phases follow immediately from 
the structure of the scattering theory. The 'true' correlation length £ is determined by the 
large distance decay of the spin-spin correlator as (aj(x)<jj(0)) ~ e - ' 11 '^. Then it follows 
from (0,0) that 

£ = £ d = l/2m (16) 

(here and below the subscript o (d) denotes quantities computed in the ordered (disor- 
dered) phase). Numerical simulations fl5| and large q expansions |16| suggest that the 



phase independence of £ holds true for all q > 4. 

Since the interfacial tension between two coexisting phases is given by the total mass of 
the lightest excitation interpolating between them, we also have a d = m and a d = cr 00 /2. 
The latter result is known to hold for all for q > 4 JIT] . 

The relation of the other interesting thermodynamic quantities (spontaneous mag- 
netisation M, latent heat L, susceptibility x, specific heat C, second moment correlation 
length ^2nd) with the connected correlators of aj and e, and their behaviour as q — > 4 + 
are 

M = (fijlo-j-m-) ~ Bf 1/8 , (17) 

L=(e) d -(e) ^Cr 1/2 , (18) 
X = J d 2 x MxfaiO)) ~ T o4 e /4 , (19) 
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d 2 x (e(ar)e(0)) ~ A o4 £ . 



(20) 



For the ordered case, the two-point correlators of Oj entering \ and ^nd are computed on 
the vacuum \ Since e(x) is odd under the duality transformation exchanging the low- 
and high-temperature phases, we have 



(e) d = -(e) a 



(21) 



The critical amplitudes are normalisation dependent but can be combined into a series 
of universal ratios characterising the scaling limit. We can evaluate the critical amplitudes 
by integrating the two-particle approximations (12 , 131) of the correlators. What we need 
to know are the two- kink form factors of the operators e and Oj. Once again the result is 



contained in Ref. 131 and reads 



with T_(0) = T_ 



2nr), 



TiL - . F (9) 



e 

4&. 



_ i e ±| (Tr+ie) 



6T.(z7r) cosh- 



T±(0)F o {0) 



(22) 
(23) 



T + (0)=exp^f- ' 



sin 



o 



F {9) 



-i sinh - exp 



x sinh 2x 

00 dx e _ 2 



x 



(2^-9)- 









sin 2 


(in — t 











(24) 
(25) 



/o x cosh | sinh a; 
The results we obtain for the universal amplitude ratios are given in Table 1 and 



compared with those following from the combination of the exact 
lattice results for the amplitudes. 
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Field theory 


Lattice 


fd 


0.6744 


0.673(8) 


foj fd 


0.9340 


0.935(5) 


ryr 


1.1406 


1.19(5) 


A d /£ 2 


0.1047 


0.105(3) 


rys 2 


0.06607 


0.0656(15) 
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Table 1. Universal amplitude ratios for the g-state Potts model at T = T c , q — > 4 + . The 
field theoretical results are obtained within the two-particle approximation. 

The accuracy exhibited by the two-particle approximation does not come as a surprise 
since it is known as a common feature of this kind of computations within integrable field 
theory. In the present case the accuracy is enhanced by the low scaling dimensions of the 
spin and energy operators which lead to mild singularities for their correlators and then 
to a small contribution of short distances to the integrals. We estimate that the errors on 
our values for the amplitude ratios do not exceed order 0.1%. 

The scaling limit we discussed so far corresponds to q — > 4 + . At q — 5, however, the 
correlation length is still some 2500 times the lattice spacing and this suggests that our 
results for q — >• 4 + could still provide the basis for an approximate description. 

For a generic value of q > 4, the model has q + 1 degenerate ground states at the 
transition point and the elementary excitations are 2q kinks going from the disordered 
vacuum to the q ordered vacua, and vice versa. If the correlation length is sufficiently large, 
an approximate scaling should still hold and then it makes sense to keep for the physical 
quantities the parameterisations (17-20), namely a power of the correlation length times 
an amplitude. It is easy to see, however, that in the present case we have to allow for a 
g-dependence of the amplitudes. Consider in fact the correlator 

G a j(x) = (Q a | (Jj (x) o-j(0)\ n«) , a = 0, i; i,j = l,...,q. (26) 

Its two-particle approximation in the disordered phase is 

G oj (x) ~ t |F ?|V^I = |i#| a e-*W , (27) 

where integration over momenta is understood. The last equality follows from colour 
symmetry and J2j a j = which imply 

F% a = {qSa - I)/ {q - l)i^. . (28) 

When integrating the correlator to obtain the amplitude of, say, the susceptibility in the 
disordered phase, the form factors computed at q = 4 should give a good approximation 
as long as £ is large. The explicit factor q/(q—l) dictated by the number of intermediate 
states and symmetry, however, has to be taken into account and is expected to determine 
the main deviation from the q — > 4 + value of the amplitude. Following the same reasoning, 
the susceptibility amplitude in the ordered phase should be basically constant in q since 
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there is only one intermediate state. More generally, we are led to expect that the ratios 
listed in Table 2 are approximately constant in q for £ large enough so that the continuum 
description is accurate. Their values determined from the results of the large q expansion 
and reported in the Table seem to confirm our picture. Our field theory results for R\ 
and R 2 as q — > 4 + are 0.855 and 0.0579, respectively. 



q 


4+ 


5 


10 


£ (lattice units) 

Jd) ,£ 

S2nd/ S 
Jo) ,Ad) 

R l = (q- l)/qXd/Xo 

r 2 = xo/my 


0.673(8) 
0.935(5) 
0.89(4) 
0.0550(6) 


2512.24 
0.671(3) 
0.934(7) 
0.810(5) 
0.0589(2) 


10.559 
0.6587(1) 
0.9579(2) 
0.80399(1) 
0.05784(1) 



Table 2. Values obtained combining the exact results of Refs. 
expansions of Ref. pi 
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Let us conclude this note by considering the 'transverse' susceptibility xt f| obtained 
integrating Gij(x) rather Gjj(x). From fl28|), in the two-particle approximation 

Xr 1 

Xo ~ (g-i) 2 ' 



(29) 



This result is basically a consequence of the nature of the elementary excitations and is 
expected to hold as a good approximation for all q > 4 at T c , as long as £ a. We are 
not aware of lattice results on xt fc> r comparison. 

In summary, we have shown that the limit q — > 4 + in the Potts model defines an 
integrable massive field theory, whose S'-matrix and form factors may be computed exactly. 
The results for integrated correlation functions are in excellent agreement with lattice- 
based numerical results. This shows how methods of continuum field theory are not 
restricted to the description of second-order transitions only. 
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